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ADDENDUM 

Propagation  in  a  Transport  System  with  60  Degree  Phase  Advance. 


The  simulations  described  in  this  report  were  performed  on  an  alternate  gradient  focusing  system 
with  a  low  current  phase  advance  of  90  degrees  per  cell,  in  this  case,  the  current  intensity  appears  to 
be  limited  by  space  charge  driven  emittance  growth.  More  recent  results  have  shown  that  a  transport 
system  with  a  60  degree  phase  advance  behaves  differently  and  the  current  does  not  appear  to  be  simi¬ 
larly  limited. 

By  considering  the  linear  stability  of  a  K-V  distribution,  Laslett  and  Smith  [1]  have  determined 
that  a  60  degree  system  can  stably  transport  a  more  intense  beam  than  a  90  degree  system.  Simulations 
similar  to  those  described  here  were  performed  to  verify  that  a  K-V  system,  with  i  60  degree  phase 
advance  depressed  to  24  degrees  by  space  charge,  will  propagate  with  no  emittance  growth  or  evidence 
of  instability. 

More  recent  simulations  performed  by  Hofmann  [2,3]  on  a  60  degree  transport  system  found  that 
a  waterbag  distribution  with  a  depressed  phase  advance  of  12  degrees  showed  no  substantial  rms  emit¬ 
tance  growth.  Motivated  by  these  results,  simulations  of  a  K-V  distribution  with  depressed  tunes  of  12 
degrees  and  6  degrees  were  run.  While  the  K-V  distribution  does  show  evidence  of  the  theoretically 
predicted  instabilities,  the  instability  saturates  at  a  very  low  level  and  rapidly  approaches  the  squared 
distribution  described  by  Hofmann.  In  contrast  to  the  90  degree  system,  very  little  rapid  growth  in  the 
rms  emittance  is  observed.  Some  very  small  slow  growth  in  the  rms  emittance  is  observed,  but  after 
propagation  through  100  magnet  pairs  the  rms  emittance  has  grown  only  about  5%.  This  absence  of 
growth  in  the  rms  emittance  is  similar  to  what  is  observed  in  the  simulation  of  continuously  focused 
channels  [3]. 

In  summary,  recent  results  have  indicated  that  beam  propagation  in  a  60  degree  channel  appears 
possible  at  much  higher  intensities  than  can  be  achieved  in  a  90  degree  system. 
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NUMERICAL  STUDY  OF  SPACE  CHARGE  LIMITED  TRANSPORT  IN  A 
QUADRUPOLE  TRANSPORT  SYSTEM  WITH  90  DEGREE  PHASE  ADVANCE 

I.  INTRODUCTION 

The  possibility  of  using  a  beam  of  energetic  heavy  ions  to  ignite  an  inertially  confined  thermonu¬ 
clear  pellet,  [1,2]  has  increased  the  interest  in  transporting  low  emittance  high  current  beams.  It  has 
been  known  for  some  time  [3]  that  such  beams  are  subject  to  space  charge  driven  instabilities.  Ques¬ 
tions  concerning  the  onset  of  these  instabilities  and  their  effect  on  beam  transport  have  been  the  sub¬ 
ject  of  a  considerable  body  [3-9]  of  theoretical  and  computational  study. 

In  this  report,  several  computer  simulations  are  presented  which  examine  the  influence  of  space 
charge  driven  instabilities  on  the  behavior  of  a  thin  lens,  alternate  gradient,  transport  system  with  a  90 
degree  phase  advance  per  cell.  These  simulations  follow  the  evolution  of  several  thousand  simulation 
particles  in  their  self  consistent  electric  fields.  By  following  the  detailed  non-linear  evolution  of  this  sys¬ 
tem  at  various  current  levels,  it  is  shown  that  these  space  charge  driven  instabilities  can  cause  emittance 
growth.  The  behavior  is  found  to  depend  on  the  current.  But  the  results,  which  show  that  increasing 
the  current  also  increases  the  emittance  growth,  suggest  that  this  emittance  growth,  by  diluting  the 
phase  space  density  of  the  beam,  may  place  an  upper  limit  on  the  current  intensity  which  can  be  tran¬ 
sported. 

In  Section  II  we  discuss  the  simulation  results,  and  in  Section  III  we  discuss  the  implications  to 
the  transport  of  intense  currents. 

II.  SIMULATION  RESULTS 

The  behavior  of  quadrupole  transport  systems  in  the  presence  of  space  charge  has  been  investi¬ 
gated  for  beams  with  Kapchinskjj- Vladimirskij  (K-V)  distributions  [11]  by  integrating  the  K-V  envelope 


Manuscript  submilled  January  13.  1981 


1 


I.  HABER 


equations  112,13],  Stability  of  K-V  systems  has  in  turn  been  investigated  by  a  perturbation  analysis  of 
these  envelope  equation  solutions  (13-17).  Regions  of  rapid  instability,  with  growth  times  of  just 

several  magnet  traversals,  have  been  discovered. 

The  envelope  equation  formalism  has  not,  however,  yielded  insight  into  the  saturation  levels  or 
the  general  non-linear  behavior  of  these  instabilities.  When  the  amplitude  of  the  unstable  modes 
becomes  large,  the  beam  can  no  longer  be  approximated  by  a  K-V  distribution.  Though  envelope  equa¬ 
tions  can  be  used  under  some  circumstances  to  predict  the  behavior  of  non  K-V  systems  [18,19],  and 
some  recent  work  has  been  done  on  instabilities  in  non  K-V  systems  [20-22],  even  the  linear  behavior 
of  these  unstable  systems  is  not  yet  understood.  Previous  simulations  [23]  have  shown  that  for 

moderate  tune  depressions,  the  system  evolves  to  a  steady  state  after  only  a  factor  of  two  growth  in  j l 

I 

emittance.  The  simulations  presented  here  were  performed  in  order  to  extend  that  work  to  a  range  of  j 

current  intensities  so  as  to  systematically  investigate  the  non-linear  evolution  of  an  unstable  K-V  sys¬ 
tem. 


The  numerical  model  is  a  two  dimensional  particle  simulation  code  of  the  particle-in-cell  type  [24] 
optimized  for  a  Texas  Instruments  ASC  computer.  The  simulation  is  performed  in  a  plane  perpendicu¬ 
lar  to  the  direction  of  beam  propagation,  in  a  reference  frame  moving  with  the  beam.  All  variation  in 
the  beam  direction  is  neglected  in  this  frame.  Since  the  transverse  behavior  is  nonrelativistic,  velocities 
and  momenta  will  be  referred  to  interchangeably  for  notational  convenience. 

Numerical  tests  have  been  performed  [2S]  which  indicated  the  existence  of  a  broad  numerical 
regime  where  the  simulation  results  are  independent  of  numerical  parameters.  These  tests  showed  that 
variations  of  such  parameters  ?s  time  step,  grid  size,  boundary  location,  and  number  of  particles  by  fac¬ 
tors  of  two  from  the  values  used,  left  the  results  essentially  unchanged.  The  simulations  described  here 
were  generally  performed  on  a  128  x  128  system  employing  16384  simulation  particles  and  using  six 
time  steps  to  follow  the  evolution  from  the  center  of  one  magnet  to  the  center  of  the  next.  In  order  to 
eliminate  the  possiblity  that  any  emittance  growth  observed  is  due  to  numerical  effects,  a  series  of 
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simulations  was  also  performed  to  test  the  dependence  of  emittance  growth  on  the  number  of  particles 
used  in  the  simulation  126]  Care  was  used  in  all  cases  to  ensure  that  the  number  of  particles  in  the 
simulations  was  sufficient  to  preclude  numerical  emittance  growth.  A  thin  lens  alternate  gradient 
transport  system  with  a  phase  advance  of  90  degrees  per  magnet  pair  has  been  assumed  for  all  the 
simulations  reported  here 

When  a  system  of  linear  lenses  is  used  to  propagate  a  beam  with  a  K-V  distribution,  the  space 
charge  forces  are  linear  and  oppose  the  net  focusing  forces.  All  the  particles  then  execute  orbits  with 
the  same  phase  advance  per  cell  and  that  phase  advance  is  reduced  from  the  low  current  value.  This 
reduced  phase  advance  is  a  measure  of  the  amount  of  current  being  transported.  In  fact,  for  a  matched 
beam  with  a  K-V  distribution  being  transported  in  a  thin  lens  alternate  gradient  transport  system,  only 
two  parameters  are  required  to  completely  specify  the  system.  So,  for  example,  specifying  the  phase 
advance  in  the  absence  of  space  charge  and  also  the  reduced  or  depressed  phase  advance  completely 
specifies  the  system.  Any  other  way  of  specifying  the  total  current  can  also  serve  instead  of  the 
depressed  tune  as  a  second  parameter. 


For  a  K-V  system  being  transported  by  finite  lenses  a  measure  of  the  transported  current  can  be 
derived  [12]  by  writing  the  envelope  equations  in  normalized  form.  For  a  thin  lens  system,  however, 
the  appropriate  [27]  measure  of  beam  intensity  is: 


Q  = 


4  q2  N  r„  f 


(1) 


A  0  y  2t 

where  q  and  A  are  the  ion  charge  state  and  mass  number,  N  is  the  number  of  particles  per  unit  length, 


rp  is  the  classical  proton  radius,  77  e  is  the  emittance,  and  /is  the  thin-lens  focal  length. 


Several  runs  have  been  performed  and  these  are  summarized  in  Table  1.  The  cases  run  are 
identified  by  the  depressed  tune,  <r.  Also  tabulated  are  the  values  of  0  as  defined  in  Eq.  1.  The  last 
column  in  Table  1  is  the  number  of  filaments.  A/,  formed  in  the  phase  space  distribution  by  the 
unstable  fields.  Note  that  the  first  two  cases  are  stable. 


3 


1  HABER 


Figure  1  is  a  plot  of  2048  sample  particles  (of  the  16K  in  the  run)  showing  four  projections  of  the 
4-dimensional  x-Px—y—Py  phase  space.  This  plot  is  of  the  initial  matched  K-V  distribution  of  one  of 
the  runs.  A  K-V  distribution  which  has  the  form 


has  a  uniform  projection  in  each  of  the  four  views  shown  in  Figure  1.  The  graininess  of  the  plots  is 
therefore  due  to  the  statistical  variations  in  the  sample. 

Figure  2  is  a  similar  plot  of  the  phase  space  of  the  45  degree  case,  but  shown  after  the  initial 
equilibrium  has  passed  through  29  magnet  pairs.  Three  well  defined  filaments  can  be  seen  in  both  the 
x—Px  and  y—Py  projections.  Phase  space  filamentation  of  this  type,  which  can  also  be  caused  by  lens 
non-linearities,  is  known  (281  to  cause  emittance  growth.  This  is  also  the  case  here.  It  is  worth  noting, 
however,  the  distinction  between  non-linearities  due  to  external  forces  and  those  seen  here.  In  the 
cases  shown  here,  the  non-linearities  grow  in  proportion  to  their  own  amplitudes  (exponential  growth), 
at  least  in  the  early  stages  of  the  instabilities.  Also  since  the  theory  of  these  instabilities,  which  exam¬ 
ines  small  perturbations  about  an  equilibrium,  does  not  predict  any  emittance  growth,  it  is  expected  that 
any  emittance  growth  will  only  occur  when  the  unstable  perturbation  is  no  longer  small.  For  a  particular 
mode  therefore,  the  observation  of  filamentation  and  emittance  growth  is  likely  to  occur  near  the  point 
where  the  growth  of  that  unstable  mode  is  saturating. 

Figure  3  is  a  plot  of  the  rms  emittances  xrmPxrm$  and  yrmsPyrms  (the  cross  terms  in  the  normal 
definition  (181  are  negligibly  small)  as  the  the  beam  evolves  from  its  initial  K-V  equilibrium.  These 
emittances  initially  remain  constant  as  the  unstable  modes  grow  from  their  small  initial  values  which 
are  non-zero  due  to  the  statistical  imperfections  in  approximating  the  initial  distribution  by  a  small 
number  (in  this  case  16384)  particles. 

Once  the  initial  perturbation  has  grown  to  a  finite  value,  the  beam  phase  space  density  rapidly 
filaments  and  the  emittance  grows  until  a  saturated  state  is  reached.  The  growth  then  stops.  The  rapid 
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growth  of  this  mode  is  in  agreement  with  the  linear  theory  of  the  stability  of  K-V  distributions  in  alter¬ 
nate  gradient  systems.  Because  only  one  unstable  mode  is  present  for  a  system  with  a  depressed  tune 
of  45  degrees,  it  has  been  possible  to  obtain  detailed  agreement  with  the  theoretical  predictions  as  to 
growth  rate  and  eigenfunction  characteristics  tor  this  case  [5i.  As  pointed  out  previously  this  agree¬ 
ment  occurs  only  before  the  major  part  of  the  emittance  growth  is  observed 

When  the  initial  current  is  increased  the  structure  of  the  most  visible  unstable  mode  changes. 
Figure  4  shows  the  phase  space  after  20  magnet  pairs  of  an  initial  K-V  distribution  with  enough  current 
to  depress  the  phase  advance  from  90  degrees  to  30  degrees  In  this  case,  four  filaments  are  visible. 
Figure  5  shows  a  continuation  of  this  trend  when  the  tune  is  depressed  further  to  15  degrees.  Finally, 
Fig.  6  shows  the  case  with  tune  depressed  to  7.5  degrees  In  this  case  the  filaments  are  not  quite  as 
well  defined  but  about  8  to  10  of  them  are  present  The  more  current  therefore,  the  greater  the 
number  of  filaments  which  are  visible.  However,  the  mode  which  causes  the  most  phase  space  struc¬ 
ture,  and  is  therefore  the  most  visible,  is  not  necessarily  the  only  one  that  is  unstable.  In  the  higher 
current  regime  it  is  in  fact  likely  that  several  modes  are  simultaneously  unstable. 

Figure  7  is  a  plot  of  the  v  and  y  rms  omittances  for  the  simulation  of  the  case  with  a  30  degree 
depressed  tune.  The  general  behavior,  i.e  .  a  period  of  constancy  followed  hy  a  period  of  rapid  growth 
and  then  a  steady  state  after  saturation  of  this  growth  is  similar  to  the  45  degree  case.  The  growth  rate 
in  this  case  is  faster.  There  is  also  an  anistropy  in  the  emittance  in  the  final  state  which  appears  to  be 
due  to  the  statistical  excitation  of  the  y-growing  eigenfunction  at  a  higher  amplitude  because  of  the  par¬ 
ticular  set  of  random  initial  conditions  on  the  particles. 

Figure  8  illustrates  the  occurrence  of  a  transition  in  behavior  when  the  tune  is  depressed  further. 
The  emittance  is  again  a  constant  initially  and  then  grows  rapidly.  However,  instead  of  reaching  a 
steady  state,  the  emittance  continues  to  grow  slowly  Since  the  distribution  function  by  this  lime  has 
evolved  far  from  the  initial  K-V  equilibrium,  it  is  difficult  to  say  whether  this  growth  is  caused  by 
unstable  modes  which  were  unstable  initially  or  modes  which  appeared  as  a  result  of  the  change  in  the 
distribution  function,  or  in  fact,  whether  the  growth  is  caused  by  the  non-linearities  of  the  distribution 
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emittance  growth,  several  numerical  tests  were  run  in  which  quantities  such  as  the  number  of  particles 
and  the  spatial  resolution  of  the  numerical  grid  were  changed.  In  none  of  these  cases  was  the  growth 
rate  of  the  emittance  affected  substantially. 

The  simulation  run  with  a  depressed  tune  of  7.5  degrees,  as  shown  in  Fig.  9,  exhibits  an  even 
more  pronounced  persistence  to  the  emittance  growth.  In  order  to  completely  rule  out  'he  possbility 
that  this  emittance  growth  was  caused  by  particle-particle  collisions,  this  case  was  run  using  t>4o00  parti¬ 
cles.  In  this  case  the  emittance  has  grown  by  about  a  factor  of  1.5  between  the  end  of  the  rapid  growth 
period  and  the  end  of  the  run.  This  greater  emittance  growth  as  the  current  increases  seems  to  indicate 
the  possibility  that  there  is  a  maximum  current  which  can  be  stably  transported.  However,  the  magni¬ 
tude  of  this  current  depends  on  the  details  of  the  scaling  of  a  beam  transport  system  and  will  be  dis¬ 
cussed  in  the  next  section. 

The  rms  emittance  is  the  product  of  a  configuration  space  and  velocity  space  component,  each  of 
which  can  evolve  separately.  Figure  10  shows,  for  the  case  with  45  degree  depressed  tune,  a  plot  of  the 
x-extent  of  the  beam  at  the  center  of  an  x  focusing  lens  as  the  system  evolves.  Three  curves  labelled  r, 
a ,  and  m  at  the  right  margin  correspond  to  the  rms,  absolute  average  and  maximum  (the  position  of  the 
maximally  displaced  particle)  measures  of  the  size. 

Figures  11-13  show  the  same  graphs  for  the  runs  with  depressed  tunes  of  30,  15.  and  7.5  degrees 
respectively  Both  the  rms  and  absolute  average  beam  sizes  change  little.  Therefore  little  of  the  growth 
in  the  rms  emittance  shown  in  Figures  3,  7,  8  and  9  occurs  due  to  growth  in  the  configuration  space 
component.  In  all  cases  the  maximum  particle,  or  edge  of  the  distribution  function,  has  moved  to 
about  twice  its  initial  value.  The  beam  size  has  not  increased  much  for  most  of  the  panicles,  but  some 
particles  have  been  scattered  out  to  twice  the  initial  edge  of  the  beam. 

Figures  14-17  are  the  corresponding  graphs  of  the  evolution  of  the  beam  extent  in  Px  space,  also 
taken  at  the  center  of  an  x-focusing  lens.  These  curves  also  illustrate  that  most  of  the  growth,  espe¬ 
cially  at  the  highest  currents,  is  due  to  growth  in  the  velocity  component  of  the  emittance.  The  edge  of 
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the  beam  has  moved  way  out  in  velocity  space.  Since  the  same  spreading  has  not  occurred  in 
configuration  space,  this  is  an  indication  of  the  highly  non-linear  space  charge  forces  near  the  edge  of 
the  beam,  with  the  non-linearity  becoming  more  pronounced  at  the  higher  currents. 

These  curves  all  refer  to  the  behavior  at  the  center  of  the  x-focusing  lenses.  This  was  done 
because  this  is  where  the  beam  is  at  its  maximum  in  the  x  direction,  and  it  is  therefore  appropriate  to 
look  at  the  beam  growth  at  this  point  in  order  to  determine  the  size  channel  which  would  be  necessary 
to  propagate  the  beam. 

The  self-consistent  equilibria  corresponding  to  a  K-V  distribution  are  not  difficult  to  visualize.  All 
projections  of  the  beam  are  uniform  ellipses.  The  non  K-V  steady  states  which  the  system  has  evolved 
to  in  these  simulations  are  somewhat  more  complex.  Figure  18  shows  four  views  of  the  x-Px—y-Py 
phase  space  from  the  simulation  with  an  initial  30  degree  depressed  tune  after  the  beam  has  propagated 
through  100  magnet  pairs.  Figures  19  and  20  are  contour  plots  of  the  x-Px  and  y—Py  phase  spaces  of 
Figure  18.  Adjacent  contours  are  separated  by  a  factor  of  four  in  density.  This  run  was  chosen  because 
it  was  the  highest  current  at  which  a  steady  state  was  reached.  The  shapes  of  the  x-Px  and  v-Py  con¬ 
tours  differ  because  the  plot  was  done  in  the  center  of  an  x-focusing  lens.  A  similar  contour  plot  done 
at  the  center  of  the  y-focusing  lens  shows  the  expected  reversal  of  the  two  plots. 

Another  set  of  plots  can  be  generated  of  the  configuration  space  and  velocity  space  contours. 
These  are  most  conveniently  taken  half  way  between  the  x-focusing  and  t-focusing  lenses.  At  this 
plane  a  matched  equilibrium  is  symmetric  in  the  x  and  y  directions.  Figures  21  and  22  show  these  plots 
for  configuration  and  velocity  space  respectively  for  the  same  case  run  with  a  30  degree  depressed  tune 
taken  at  the  end  of  the  simulation. 

The  squaring  of  these  contour  lines  is  just  one  indication  as  to  the  difficulty  of  analytic  characteri¬ 
zation  of  the  high  space  charge  equilibria.  This  squaring  increases  as  the  current  is  increased.  Since  the 
contours  in  phase  space  are  not  ellipses,  it  is  also  an  indication  that  any  definition  of  emittan ce  in  terms 
of  the  percentage  of  particles  within  an  elliptical  phase  space  contour  may  not  be  appropriate. 
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For  completeness,  it  is  also  possible  to  plot  various  cross  sections  of  these  curves.  Figure  23  is 
such  a  cross  section  of  Figure  21.  Figure  23  is  a  plot  of  fix),  the  number  density  as  a  function  of  x,  at 
the  y  =  0  plane  integrated  over  all  velocities.  Because  this  curve  is  generated  at  the  symmetry  plane,  it 

is  the  same  form  as  a  curve  of  f(y)  at  x  -  0.  Figure  24  is  a  similar  plot  of  f(Py)  at  Px  =  0  integrated 
over  all  positions.  These  curves  are  typical.  A  reasonable  approximation  to  them  can  be  obtained  as 

the  sum  of  a  distribution  with  a  parabolic  profile  and  a  smaller  one  with  a  linearly  decreasing  tail  (21] 
The  number  of  such  curves  which  can  be  generated  for  all  the  cases  run  is  large.  Most  of  the  informa¬ 
tion  contained  in  these  curves  can  be  summarized  by  noting  that  as  the  current  increases,  the  contours 
get  squared  as  do  the  configuration  space  cross  sections.  The  velocity  space  cross  sections  develop 
increasingly  long  tails.  The  information  of  the  length  of  these  tails  however,  can  be  pretty  satisfactorily 
estimated  from  the  curves  of  the  evolution  of  velocity  vs.  time  (Figs.  14-17),  by  comparing  the  max¬ 
imum  velocity  particle  to  the  rms  velocities.  However  at  the  mid  plane  the  velocity  and  configuration 
behavior  look  generally  similar  to  Figs.  23  and  24. 

III.  CONCLUSIONS 


Current  intensity  depends  on  the  ratio  of  current  to  emittance.  Therefore  as  the  current 
increases,  if  the  emittance  increases  correspondingly,  then  the  depressed  tune  is  unchanged.  In  fact,  if 
the  distribution  function  remains  unchanged,  the  beam  is  then  represented  by  the  same  solution  to 
Maxwell’s  equations  in  the  transverse  direction.  It  is  possible  from  purely  dimensional  arguments  to 
show  that  any  solution  to  Maxwell’s  equations  in  the  transverse  direction  corresponds  to  a  famiK'  of 
solutions  at  different  currents.  These  currents  can  be  related  to  the  power  propagating  in  a  transport 
system  by  the  relation 


-  C  — 
z 


4/3 

(Bq  t)i,}  (fiy)Vi  (y- 1)  (3) 

where  P  is  the  power  in  watts,  A  and  z  are  the  mass  number  and  charge  stale.  Bq  is  the  quadrupole  pole 
tip  field  in  Tessla,  and  rr  <  is  the  emittance  in  meter-radians.  The  coefficient  C  contains  the  informa¬ 
tion  on  the  current  intensity  and  the  beam  distribution  function. 
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This  formula  was  first  suggested  by  Maschke  and  was  derived  by  Courant  [29]  for  a  K-V  distribu¬ 
tion.  A  different  treatment  was  derived  by  Reiser  [30].  The  derivation  in  terms  of  purely  dimensional 
considerations  is  presented  in  the  appendix. 

It  is  possible  to  determine  a  value  for  the  coefficient  C  in  Eq.  (3)  corresponding  to  any  solution  to 
Maxwell’s  equations  in  general,  and  therefore  specifically  for  the  distributions  coresponding  to  each  of 
the  simulations  performed.  However,  the  power  transported  is  expressed  in  Eq.  (3)  in  terms  of  the 
pole  tip  field.  Therefore  because  the  simulations  were  performed  using  thin  lenses,  and  these  are 
unphysical  as  they  require  singular  pole  tip  strengths,  some  assumption  must  be  made  about  the  magnet 
fill  in  the  transport  system.  From  the  extensive  linear  stability  calculations  [13-15],  it  is  known  that,  at 
least  in  the  linear  stages  of  the  instability,  the  behavior  of  the  system  is  largely  insensitive  to  magnet  fill 
as  long  as  the  phase  advance  per  cell  is  kept  unchanged.  This  insensitivity  to  lens  fill  has  also  been 
seen  in  simulations  [7], 

Laslett  has  suggested  [31]  the  definition  of  a  beam  figure  of  merit,  /  =  (?/<wmax)2/\  based  on  the 
behavior  of  the  scaled  envelope  equation  [12].  The  coefficient  C  in  Eq.  (3)  is  related  to  F  by,  C  =  3.43 
x  1015  fin  MKSA  units.  For  a  K-V  system,  the  values  of  Fhave  been  calculated  for  a  range  of  param¬ 
eters  [31-33].  These  results  are  plotted  in  Fig.  25,  as  a  function  of  depressed  tune,  for  a  system  with 
90  degree  phase  advance  and  with  50  percent  lens  fill.  The  values  of  /  can  be  calculated  beyond  the 
range  shown  in  Fig.  25  from  the  expression  [31]  F  =  9.5  o-_2/3,  which  is  accurate  to  within  one  percent 
below  10  degree  depressed  phase  advance.  From  the  K-V  values,  it  is  possible  to  calculate  the  value  of 
F,  and  therefore  the  coefficient  C,  appropriate  to  the  final  state  in  each  of  the  simulation  runs. 

Since  the  current  does  not  change  during  the  simulation,  the  coefficient  in  Eq.  (4)  changes  due  to 
the  growth  in  emittance  as  the  beam  evolves.  Because  the  beam  size  also  increases  an  additional  factor 
must  be  included  to  account  for  the  fact  that  the  magnet  pole  lip  must  be  moved  out.  These  values  are 
listed  in  Table  II.  The  value  at  60  degrees,  which  is  not  unstable,  is  the  same  as  Laslett's  value.  The 
other  values  are  calculated  from  Laslett's  values  by  dividing  by  the  increase  in  rms  emittance  to  the  2/3 
power  and  then  also  dividing  by  the  increase  in  size  to  the  beam  edge  (as  measured  from  Fig.  13-16) 
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also  raised  to  the  2/3  power.  This  second  factor  is  necessary  if  the  pole  tip  is  moved  out  far  enough  to 
include  the  entire  beam.  Since  the  emittances  used  are  rms  cmittance  these  particular  coefficients  are 
valid  only  when  using  this  rms  measure. 

The  last  two  entries  in  the  table  correspond  to  distributions  which  are  not  in  steady  state,  but  are 
still  evolving.  Notice  for  example  that  the  initial  figure  of  merit  increases  by  a  factor  of  1.6  between  15 
and  7.5  degrees,  by  the  end  of  the  run  this  improvement  has  been  reduced  to  1,2.  If  the  run  were  con¬ 
tinued  further,  this  reduction  would  have  continued  still  further  resulting  in  the  possibility  of  no  gain  in 
transported  current  intensity. 

For  the  simulations  performed,  that  is  a  constant  current  propagating  down  a  uniform  transport 
line,  a  K-V  system  can  become  rapidly  unstable.  This  rapid  instability  appears  to  saturate  after  about  a 
factor  of  two  growth  in  cmittance.  At  the  highest  current  intensities  studied,  this  period  of  rapid 
growth  is  followed  by  a  period  of  slower  growth  which  tends  to  limit  the  current  intensity  than  can  be 
propagated  over  long  distances. 

Because  of  the  inherently  non-linear  nature  of  the  phenomena  involved  it  is  difficult  to  say  any¬ 
thing  about  beam  distributions  more  general  than  K-V  If  it  as  assumed  that  the  behavior  of  non  K-V 
distributions  is  adequately  represented  by  the  behavior  of  the  simulations  alter  initial  saturation,  then  it 
would  be  possible  to  propagate  very  intense  currents  for  modest  distances  without  too  much  emiilance 
growth  However  the  existence  of  modes  which  can  grow  rapidly  and  which  could  be  excited  if  the 
beam  distribution  would  evolve  to  a  state  closer  to  a  K-V  distribution,  suggests  that  caution  be  exer¬ 
cised.  Conservative  design  still  seems  to  preclude  operation  in  this  space  charge  limited  regime  without 
a  much  fuller  understanding  of  the  effect  of  beam  distribution  function  on  stability. 
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Table  I.  Summary  of  Cases  Run 
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Table  II.  Transport  Figures  of  Merii  Modified 
by  Space  Charge  Driven  Instabilities 
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Appendix 

CURRENT  SCALING  IN  A  STRONG  FOCUSING  CHANNEL 


Consider  a  beam  of  charged  particles,  in  the  paraxial  approximation,  in  a  reference  frame  moving 
with  the  beam.  The  transverse  dynamics,  in  this  frame,  is  described  by  the  Vlasov-Poisson  system  of 
equations. 

|^  + v-V/+ (E„  +  E)  •  Vv/=  0  (Al) 

where  the  distribution  function  is  normalized  so  that 

J*  /  </v,  dvv  dx  dy  =  1 . 

and  E,,v  represents  the  focusing  magnet  force  in  the  beam  reference  frame.  The  space  charge  fields  are 
found  from 


V  •  E  »*  4 irp. 


( A2) 


In  order  to  investigate  the  scaling  of  this  system  we  write  these  equations  in  terms  of  normalized 
variables.  Setting  r  =  r(lr,  v  =  v„  v  and  t  =  t„t  with  r„  =  r0 ,  and  also  E,,A  =  E0  E„, 


H  +  v  •  V/+  ~  (£„  E„,  +  E)  •  V,  /  =  0. 
dr  v.fm 


(A3) 


Since  the  divergence  operator  is  just  a  linear  operator  for  the  purpose  of  this  scaling  analyses,  it 
can  be  formally  inverted  and  written  as: 


E  =  4t rLp  (A4) 

where  L  is  just  the  appropriate  linear  operator  inverse  to  the  divergence.  The  self-consistent  electric 

field  can  then  be  written  in  terms  of  the  total  current  I  as 


E  = 


47T  /  r„  L 

rjfiyc 


f  f  d\\  d\\  . 


( A5) 


16 


NRL  MEMORANDUM  REPORT  4474 


Equation  (A3)  may  then  be  rewritten  as: 


U  +  V -V/  +  (C,  £„,  +  C'2£  /  /  </vv  dvy)  ■  7,  /  =  0 


where 


c,  -  -r-  £„ 


Pycvjm 


Any  distribution  function  /satisfy  (A6)  with  a  particular  form  of  £(,v,  will  therefore  be  a  solution 
for  any  combination  of  parameters  which  leave  C'i  and  C2  unchanged.  We  can  use  this  to  derive  a  scal¬ 
ing  relation  for  /. 


Introducing  an  emittance 


«  =  C’j  4’„/)3yT, 


/nCi  i'„3  w 


Solving  for  the  current 


—  E,P  f3i/}  yV3 


Cf ’  4tr  U 


If  the  external  focusing  force  is  written  in  terms  of  transverse  magnetic  fields,  E„  =  /3yfi0,  this 


becomes 


/  =  <\|^  ’C 2/}  B„vi 


The  power  transported  by  this  current,  is 


£  =  (y—  1 )  —  / 
<1 


6 2/3  fl/3  (0y)7/3  (y—  I ) 
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Fig  >  —  Four  projections  of  the  four-dimensional  K-V  distribution  showing  2048  sample  particles.  The  configuration 
space  units  are  numerical  grid  cell  and  the  velocity  (momentum)  units  are  in  grid  cells  traversed  in  one  magnet  period. 
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I  ig  1  -  Evolution  of  the  x  and  v  rms  emlttances  in  a  system  with  an  initial  matched  K-V  distribution  with  enough 
current  to  depress  the  tune  from  90  degrees  to  45  degrees.  The  vertical  scale  is  normalized  to  the  initial  emittance. 
The  horizontal  scale  is  measured  in  alternate  gradient  magnet  pairs. 


f  ig  -I  —  Four  pha.se  space  projections  showing  the  evolution,  after  20  magnet  pairs,  of  a  K-V  system 
with  a  depressed  phase  advance  of  30  degrees 
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Evolution  of  the  x  and  y  rms  emittances  in  a  system  with  an  initial  K-V  distribution  with  a  depressed  phase 

advance  of  30  degrees 


Evolution  of  the  x  and  y  rms  emittances  in  a  system  with  an  initial  depressed  phase  advance  of  15  degrees. 
Some  slow  long  time  scale  growth  in  the  emittance  is  apparent. 


Evolution  of  the  x  and  y  rms  emittances  in  a  system  with  an  initial  depressed 
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Evolution  of  the  beam  size  in  the  x  direction  for  an  initial  15  degree  phase  advance 


simulation  vuth  a  4^  degree  initial  phase  advance 


f  volution  of  the  v-momenta  for  the  case  with  ^  5  degree  phase  advance 


Ill  -*« 


Fig  in  -  Four  views  of  the  x-Px-y—Py  phase  space  of  a  simulation  with  an  initial  depressed  phase  advance  of  30  de¬ 
grees  which  has  been  propagated  through  100  magnet  pairs.  The  system  has  gone  unstable  and  by  this  point  has 
reached  a  steady  state 
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A  contour  plot  of  the  v-  l‘\  phase  space  projection  shown  in  f  ig  18 

factor  of  four  in  density. 


Adjacent  countours  are  separated  by  a 
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A  conlour  plot  of  (he  density  in  real  space  of  the  simulation  with  30  degree  depressed  phase  advance  after 
100  magnet  pairs  The  plot  is  taken  midway  between  an  r-locusing  and  r-focusing  lens 
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